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Abstract. Registration, which aims to find an optimal 1-1 correspondence between shapes, is 
an important process in different research areas. Conformal mappings have been widely used to 
obtain a diffcomorphism between shapes that minimizes angular distortion. Conformal registrations 
are beneficial since it preserves the local geometry well. However, when landmark constraints are 
enforced, conformal mappings generally do not exist. This motivates us to look for a unique landmark 
matching quasi-conformal registration, which minimizes the conformality distortion. Under suitable 
condition on the landmark constraints, a unique diffcomporphism, called the Teichmiiller extremal 
mapping between two surfaces can be obtained, which minimizes the maximal conformality distortion. 
In this paper, wc propose an efficient iterative algorithm, called the Quasi-conformal ( QC) iterations, 
to compute the Teichmiiller mapping. The basic idea is to represent the set of diffcomorphisms 
using Beltrami coefficients (BCs), and look for an optimal BC associated to the desired Teichmiiller 
mapping. The associated diffcomorphism can be efficiently reconstructed from the optimal BC 
using the Linear Beltrami Solver(LBS). Using BCs to represent diffcomorphisms guarantees the 
diffcomorphic property of the registration. Using our proposed method, the Teichmiiller mapping can 
be accurately and efficiently computed within 10 seconds. The obtained registration is guaranteed 
to be bijective. The proposed algorithm can also be extended to compute Teichmiiller mapping 
with soft landmark constraints. We applied the proposed algorithm to real applications, such as 
brain landmark matching registration, constrained texture mapping and human face registration. 
Experimental results shows that our method is both effective and efficient in computing a non- 
overlap landmark matching registration with least amount of conformality distortion. 

Key words. Teichmiiller extremal mapping, quasi-conformal mapping, Beltrami coefficient, 
Linear Beltrami Solver, Landmark matching registration 

1. Introduction. Registration refers to the process of finding an optimal one- 
to-one correspondence between images or surfaces. It has been extensively applied 
to different areas such as medical imaging, computer graphics and computer visions. 
For example, in medical imaging, registration is always needed for statistical shape 
analysis, morphometry and processing of signals on brain surfaces (e.g., denoising 
or filtering). While in computer graphics, surface registration is needed for texture 
mapping, which align each vertex to a position of the texture image, to improve the 
visualization of the surface mesh. Developing an effective algorithm for registration 
is therefore very important. 

Conformal mappings have been widely used to obtain smooth 1-1 correspondences 
between different domains that minimize angular distortions. Conformal mappings 
are advantageous since it preserves the local geometry well. According to the con- 
formal geometry, given two simply-connected domains, there always exists a unique 
conformal mapping between them up to a Mobius transformation. However, the exis- 
tence of conformal mappings cannot be guaranteed in general situations. For example, 
conformal mapping between two multiply-connected domains usually does not exist. 
Besides, in practical situation, obtaining a registration that matches landmark fea- 
tures consistently arc often required. When landmark constraints are enforced, the 
existence of conformal mappings cannot be guaranteed. This motivates us to look for a 
unique landmark-matching registration, which minimizes the conformality distortion 
as much as possible. 

Under suitable conditions on the landmark constraints, there exists a unique 
quasi-conformal mapping, called the Teichmiiller extremal mapping, which minimizes 
the maximal conformality distortion. To compute this Teichmiiller mapping, we pro- 
pose in this paper an efficient and effective iterative algorithm, which is called the 
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Quasi-conformal (QC) iterations. The basic idea is to represent the set of diffeomor- 
phisms using Beltrami coefficients (BCs), and look for an optimal BC associated to 
the desired Teichmuller mapping. The associated diffeomorphism can be efficiently 
reconstructed from the optimal BC using the Linear Beltrami Solver(LBS). Given a 
set of landmark constraints, the algorithm is able to determine an optimal 1-1 cor- 
respondence (including the boundary correspondence in the case of open surfaces) 
between shapes automatically, which minimizes the conformality distortion. Besides, 
the proposed algorithm can also compute extremal mapping with soft landmark con- 
straints. It becomes necessary when landmark features cannot be accurately located, 
and hence it is better to compute registration with landmarks approximately (but 
not exactly) matched. Another major advantage of using Teichmuller mappings for 
landmark matching registrations is that the bijectity (1-1, onto) of the registrations 
can be guaranteed. Obtaining a bijective landmark matching registration is gener- 
ally difficult, especially when a large number of landmark constraints are enforced. 
Using our proposed method, a bijective Teichmuller mapping can be accurately and 
efficiently computed within 10 seconds. To test the effectiveness of our method, we 
applied the proposed algorithm to real applications, such as brain registration, con- 
strained texture mapping and human face registration. Experimental results shows 
that our method is both effective and efficient in computing a non-overlap landmark 
matching registration with least conformality distortion. 

In short, the contributions of this paper are three-folded. Firstly, we propose an 
efficient algorithm for obtaining the unique Teichmuller extremal mapping between 
shapes with landmark constraints enforced. The mapping is guaranteed to be bi- 
jective and minimizes the maximal conformality distortion. Secondly, we propose 
an algorithm to compute Teichmuller extremal mapping with soft landmark con- 
straints. Landmarks are not exactly matched, but less conformality distortion will be 
introduced. Thirdly, we apply the proposed algorithms to real applications, namely, 
constrained texture mapping, medical image registration and human face registration. 

2. Previous work. In this section, we will introduce some previous works closely 
related to our paper. 

Surface parameterization and registration have been extensively studied, for which 
different kinds of bijective surface maps have been proposed. Conformal registration, 
which minimizes angular distortion, have been widely used to obtain a smooth 1-1 
correspondence between surfaces O [3 El El [TTJ [8] . For example, Hurdal et al. [11] pro- 
posed to compute the conformal parameterizations using circle packing and applied it 
to registration of human brains. Gu et al. [3[9, 8 a proposed to compute the conformal 
parameterizations of Riemann surfaces for registration using harmonic energy mini- 
mization and holomorphic 1-forms. Later, the authors proposed the curvature flow 
method to compute conformal parameterizations of high-genus surfaces onto their 
universal covering spaces |32, 3S1 |34j . Conformal registration is advantageous for it 
preserves the local geometry well. 

Sometimes, deformations between objects might not not conformal. Instead, cer- 
tain amount of angular distortion could be introduced. To tackle with this situation, 
quasi-conformal mappings are proposed to obtain smooth 1-1 correspondence with 
bounded conformality distortion. Lui et al. [18] proposes to compute quasi-conformal 
registration between hippocampal surfaces which matches geometric quantities (such 
as curvatures) as much as possible. A method called the Beltrami Holomorphic flow 
is used to obtain the optimal Beltrami coefficient associated to the registration [T7] . 
Beltrami coefficient has been applied to represent general surface homeomorphisms, 
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which is comparatively easier to manipulate than 3D coordinate functions. Using 
Beltrami representation, compression of surface maps has been proposed [15], which 
can be applied for video compression [5D]. Wei et al. [3T] also proposes to compute 
quasi-conformal mapping for feature matching face registration. The Beltrami coef- 
ficient associated to a landmark points matching parameterization is approximated. 
However, either exact landmark matching or the bijectivity of the mapping cannot 
be guaranteed, especially when very large deformations occur. In order to compute 
quasi-conformal mapping from the Beltrami coefficients effectively. Quasi- Yamabe 
method is introduced, which applies the curvature flow method to compute the quasi- 
conformal mapping H5j. The algorithm can deal with surfaces with general topolo- 
gies. Later, quasi-conformal mapping, which minimizes the conformality distortion, 
has been studied. Zorin et al. [36 proposes an algorithm to compute the extremal 
quasi-conformal mapping between connected domains with given Dirichlet condition 
defined on the whole boundaries. The extremal mapping is obtained by minimizing 
a least square Beltrami energy, which is non-convex. The algorithm can obtain an 
extremal mapping when initialization is carefully chosen. However, the convergence 
to the global minimum cannot be guaranteed. 

Most of the above registration algorithms using conformal and quasi-conformal 
mappings cannot match feature landmarks, such as sulcal landmarks on the human 
brains, consistently. To alleviate this issue, landmark-matching registration algo- 
rithms are proposed by various research groups. Wang et al. [12 E3 E2 HH EH ES] 
proposed to compute the optimized conformal parameterizations of brain surfaces by 
minimizing a compounded energy [TSJ [15] . The obtained registration can obtain an 
optimized conformal map that better aligns the features, however, landmarks cannot 
be exactly matched. Besides, bijectivity cannot be ensured when large number of 
landmark constraints are enforced. To solve this problem, smooth vector field has 
also been applied to obtain surface registration. Lui et al. (T3l [14] proposed the use 
of vector fields to represent surface maps and reconstruct them through integral flow 
equations. They obtained shape-based landmark matching harmonic maps by look- 
ing for the best vector fields minimizing a shape energy. The use of vector fields to 
compute the registration makes optimization easier, although it cannot describe all 
surface maps. An advantage of this method is that exact landmark matching can be 
guaranteed. Time dependent vector fields can also be used. For example, Joshi et 
al. |10j proposed the generation of large deformation diffeomorphisms for landmark 
point matching, where the registrations are generated as solutions to the transport 
equation of time dependent vector fields. The time dependent vector fields facilitate 
the optimization procedure, although it may not be a good representation of surface 
maps since it requires more memory. The computational cost of the algorithm is also 
expensive. 

3. Mathematical Background. In this section, we describe some basic math- 
ematical concepts related to our algorithms. For details, we refer the readers to 

A surface S with a conformal structure is called a Riemann surface. Given two 
Ricmann surfaces M and N, a map / : M — > N is conformal if it preserves the 
surface metric up to a multiplicative factor called the conformal factor. An immediate 
consequence is that every conformal map preserves angles. With the angle-preserving 
property, a conformal map effectively preserves the local geometry of the surface 
structure. 

A generalization of conformal maps is the quasi-conformal maps, which are ori- 
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Fig. 2.1. (A) shows a human face with circle packing texture. Under the conformal param- 
eterization, infinitesimal circles are mapped to circles as shown in (B). Under quasi- conformal 
parameterization, infinitesimal circles are mapped to ellipse as shown in (C). (D) illustrates how 
the Beltrami coefficient measure the conformality distortion of a quasi- conformal map. 



cntation preserving homeomorphisms between Riemann surfaces with bounded con- 
formality distortion, in the sense that their first order approximations takes small 
circles to small ellipses of bounded eccentricity [S]. Mathematically, /: C — > C is 
quasi-conformal provided that it satisfies the Beltrami equation: 

oz oz 

for some complex valued function ft satisfying ||/i||oo < 1- A 1 is called the Beltrami 
coefficient, which is a measure of non-conformality. In particular, the map / is con- 
formal around a small neighborhood of p when ft(p) = 0. Infinitesimally, around a 
point p, f may be expressed with respect to its local parameter as follows: 

f(z) = f(p) + f z {p)z + M P )z 

= /(p)+/»(p)(* + Mp)Z). 

Obviously, / is not conformal if and only if fi(p) ^ 0. Inside the local parameter 
domain, / may be considered as a map composed of a translation to f(p) together 
with a stretch map S(z) = z + fJ,(p)z, which is postcomposed by a multiplication of 
f z (p), which is conformal. All the conformal distortion of S(z) is caused by fJ>(p). 
S(z) is the map that causes / to map a small circle to a small ellipse. From fj,(p), 
we can determine the angles of the directions of maximal magnification and shrinking 
and the amount of them as well. Specifically, the angle of maximal magnification 
is arg(/i(p))/2 with magnifying factor 1 + |/i(p)|; The angle of maximal shrinking is 
the orthogonal angle (arg(/i(p)) — ir)/2 with shrinking factor 1 — |/x(p)|. Thus, the 
Beltrami coefficient ft gives us all the information about the properties of the map 
(See Figure gD)). 

The maximal dilation of <f> is given by: 

i + IMU (33) 

1 - IMloo 



Let f = u + y/—lv. From the Beltrami equation (3.1 ), 



U) = (u x - Vy ) + ^Uv x + u v ) 

[U X + Vy) + V-1(V X - Uy) 
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<p a (U a ) C <h {Up) 

Fig. 3.1. Illustration of how Beltrami differential is defined on general Riemann surfaces. 




Fig. 3.2. Illustration of quasi-conformal mapping between Riemann surfaces. 



Let p(f) = p + \f^-i t. We can write v x and v y as linear combinations of u x and 



-v y = otiu x + a 2 u y ; 
v x = a 2 u x + a 3 u y . 



(p-_ii 2 +rf 

Similarly, 



—u y = ol\V x + a 2 v y ; 
u x = ca 2 v x + a 3 v y . 



Since V • [ y ) = 0, we obtain 

v x 



V- [A 



and V • [A 



= 



(3.5) 



(3.6) 



(3.7) 
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Quasiconformal mapping between two Riemann surfaces R\ and R2 can also be 
denned. Instead of the Beltrami coefficient, the Beltrami differential has to be used. 
A Beltrami differential fJ-(z)^ on the Riemann surface Ri is an assignment to each 
chart (U a ,4> a ) of an L^ complex- valued function fj, a , defined on local parameter z a 
such that 

on the domain which is also covered by another chart {Up, (ftp), where = -^-(f> a /3 



and 4> a fs — 4>p o <f) a l (See Figure 3.1) 



An orientation preserving diffeomorphism / : R\ — > R2 is called quasi-conformal 

d: 
dz 

, r, £ r, f 1 to mi aei_nAyifArTnol acoApioforl niif In m ( *y ^ 

dz a 



associated with /i{z)^ if for any chart (U a , <p a ) on i?i and any chart (Vp, ipp) on R2, 
the mapping f a p := ipp o / o f~ x is quasi-conformal associated with ^a{z a )^- (See 



Figure [3T2| 

Now, suppose i?i and R2 are open Riemann surfaces with the same topology. The 
boundary dilation of <p is the infimum of the maximal dilation of ?p\u over all 

quasi-conformal maps ip isotopic to cf> relative to the boundary and all neighborhoods 
U of dRt. 

In case Ri and R2 are simply-connected, conformal mapping between i?i and R2 
always exists. However, conformal mapping may not exist between multiply-connected 
domains. For example, there is generally no conformal mapping between two annuli 
with different radii of inner circles. One might be interested in investigating extremal 
quasiconformal mappings, which is extremal in the sense of minimizing the 1 1 • | |oo over 
all Beltrami differentials corresponding to quasi-conformal mappings between R\ and 
R 2 . Extremal mapping always exists but need not be unique. More specifically, an 
extremal quasi-conformal mapping can be defined mathematically as follows: 

Definition 3.1. Let f : R\ — > R2 be a quasi-conformal mapping between R\ and Rq. 
f is said to be an extremal mapping if for any quasi-conformal mapping h : R± — > R2 
isotopic to <fi relative to the boundary, 

K{4>) < K(<ip) (3.9) 



It is uniquely extremal if the inequality (3.9) is strict 



Another kind of mapping, called the Teichmuller mapping, is closely related to 
extremal mapping. Teichmuller mapping is defined as follows: 



Definition 3.2. Let f : R± — > R2 be a quasi-conformal mapping, f is said to be a 
Teichmuller mapping associated with ip : R\ — > C if its associated Beltrami coefficient 
is of the form: 

/*(/) = ^0 ( 3 - 10 ) 

for some constant k < 1 and ip =/= 0. 

Extremal mapping is not unique for general cases. However, a Teichmuller map- 
ping associated with a holomorphic function is the unique extremal mapping in its ho- 
motopic class. The Strebel's theorem explains the relationship bewtween Teichmuller 
mapping and extremal mapping. 
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Fig. 3.3. Difference between a general QC map and a Teichmiiller map. (A) shows the original 
textured mesh. It is mapped to another disk by a general QC map. Note that the distribution of the 
norm of BC are spread out. (C) shows the Teichmiiller map, whose BC norm is concentrated near 

04- 



Theorem 3.3 (Strebel's theorem). Let f be an extremal quasi- conformal map- 
ping with K(f) > 1. If Ki[f] < K{f), then f is a Teichmuller map associated with an 
integrable holomorphic quadratic function on Ri . Hence, f is also an unique extremal 
mapping. 

In particular, Teichmuller mapping and extremal mapping of the unit disk are 
closely related. 

Theorem 3.4. Let g : <9B — > <9B be an orientation-preserving homeomorphism 
of dD. Suppose further that h'(e te ) ^ and h" ie %e ) is bounded. Then there is a 
Teichmuller mapping f of finite norm that is the uniquely extremal extension of g to 
D. That is, f : D — » D is an extremal mapping with f\dn = .9- 

In other words, an extremal mapping of the unit disk with suitable boundary 
condition is a Teichmuller mapping. It can thus be obtained by searching for an 
optimal Beltrami coefficient whose maximal dilatation is the minimum, while its norm 
is constant everywhere. It turns out that in most situations, extremal quasi-conformal 
mapping is a Teichmuller mapping (even for domains with non-trivial topologies) . In 
some rare situations when an extremal mapping is not exactly a Teichmuller mapping, 
one can get a Teichmuller mapping whose dilation is arbitrarily close to the extremal 
dilation. 

Theorem 3.5. Let $ be a class of quasi-conformal mappings between the open 
Riemann surfaces R\ and R2, which are homotopic modulo the boundary. Let Kq 
be the smallest maximal dilation of the mappings in ft. Then there are Teichmuller 
mappings in associated with a meromorphic function with at most one simple pole, 
whose dilation is arbitrarily close to Kq. 

4. Mathematical Formulation. In this section, we give the mathematical for- 
mulation for obtaining the Teichmiiller extremal mapping with least conformality 
distortion. We propose to use the Beltrami coefficient(BC) to represent the mapping, 
instead of the commonly used representations by deformation fields or coordinate 
functions. The difieomorphic property of the registration can then be effectively con- 
trolled. Our goal is to formulate the problem into a variational problem to obtain an 
optimal BC, /i(/), associated to the desired extremal mapping /. 
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Suppose Di and D 2 are two domains in the complex plane with the same topol- 
ogy. Di and D 2 can either be simply-connected or multiply-connected. Suppose the 
boundary condition of the desired extremal mapping / : D\ — > D 2 is known. De- 
note it by f\dD! '■ dDi — > dD 2 = g. The Teichmiiller extremal mapping can be 
mathematically described as follows: 

% = k ri^~ and f\d D i = 9 on dD i (4- 1 ) 
oz \ip\ dz 

for some constant k and holomorphic function ip : D\ — > C. 

Recall that a Teichmiiller extremal mapping is extremal in the sense of minimizing 
the || • I |oo over all Beltrami differentials corresponding to quasiconformal mappings 
in the Teichmiiller equaivalence class. In other words, for any h : D\ — > D2 satisfying 
^■IdDi = 9, we have 

IKOIloo < IMJOHoo (4.2) 

where fJ.(f) and ^(h) are the Beltrami coefficient of / and h respectively. Hence, our 
original problem (|4.1| can be formulated as a variational problem as follows: 



/ = argmin /:£)l _ ! . r , 2 £a(/) 
:= argmin^^JIK^IU + ||V |/x(/)| |||} 

subject to: 

• fldDt — 9 (boundary condition); 

• m(/) = f° r some constant k and holomorphic function ip : Di — > C. 

Theoretically, a diffeomorphism / is associated to a unique smooth BC /i(/) with 
IIm(/)IIoo < 1- The Beltrami coefficient 11(f) measures the conformality distortion 
of the map /. It can be considered as a unique representation of /. The first term 
of the energy functional Ei aims to minimize the maximal conformality distortion 
of the mapping. The second term minimizes the harmonic energy of |/i(/)|. Since a 
Teichmiiller map has a constant norm, minimizing the second term aims to obtain a 
optimal BC whose norm is constant. 

However, minimizing E\(f) with respect to the space of diffeomorphisms between 
Di and D 2 is difficult. 



/ = argmin / {||/i(/)|| 0O + ||V K/)| 



df/dz.. ,df/dz. 2 (4.4) 



argmi n/ ||-— -|U + ||V 



df/dz df/dz 



subject to f\dD t — g and 11(f) = k— for some constant k and holomorphic function 
if : D x C. 

In order to minimize the above constrained minimization effectively, we propose 
to reformulate the energy functional with respect to space of all Beltrami coefficients: 

{v, f) = argmin^^c^M ^ g . 

:= argmin I7!iJl ^. c {||i/|| 00 + ||V \v\ || 2 } 
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subject to: 

• v = n(J) and |M|oo < 1; 

• v = k^ for some constant k and holomorphic function ip : D\ — > C; 



f\dDi =9 (boundary condition). 



In other words, the minimization problem (4.3) is reformulated to be optimized 
with respect to BCs, which are complex- valued functions defined on Di. Minimizing 
the energy functional with respect to BCs is advantageous since the diffcomorphic 
property of the mapping can be easily controlled. Every diffeomorphism is associated 
to a smooth Beltrami coefficient n(f). £*(/) measures the bijectivity (1-1 and onto) 
of /. In fact, fJ,(f) is related to the Jacobian J(f) of / by the following formula: 

W)| 2 = l^| 2 (l-K/)| 2 ) (4-6) 
oz 

Therefore, the map / is bijective if |/x(/)| is everywhere less than 1. When solving 



the minimization problem (4.5), the bijectivity of the mapping in each iterations can 
be ensured by enforcing ||^||oo < 1- Our goal is to look for an optimized smooth BC, 
v, such that its associated quasi- conf or mal map is our desired Teichmiiller extremal 
mapping. 



The boundary condition in the variational problem (4.5) can be relaxed. The 
Dirichlet condition defined on the whole boundary is not required. Also, interior 
landmark constraints can be enforced. Our goal is to solve the variational problem 
with these landmark constraints, which determine the optimal 1-1 correspondence 
(including the boundary correspondence) automatically. In other words, the boundary 



condition in the problem (4.5) can be reformulated as 



f(ai) = bi; f{pj)=q 3 ; for % = 1, ...,n; j = 1, ...,m (4.7) 

where di and hi are corresponding landmark points or curves defined on dD\ and 
dD 2 respectively; and pj and qj are corresponding interior landmark points of curves 
in Di and D 2 respectively. By optimizing the energy functional |4.5| an optimized 
Teichmiiller extremal mapping can be obtained, which matches landmark features 
consistently while minimizing the maximal conformality distortion. 

Note also that the above formulation is designed for extremal mapping of 2D con- 
nected domains. However, it can easily be extended to simply-connected or multiply- 
connected open surfaces. Let Si and S 2 be two connected open surfaces with the 
same topology. We can conformally parameterize Si and S 2 by </>i : Si —> D\ C C 
and <f>2 '■ S2 D 2 C C respectively. Then the extremal mapping / : Si S 2 between 
Si and 1S2 induces an extremal mapping / := 4> 2 o / o (jy^ 1 : D\ — > D 2 . All the above 
formulation applies to /. In other words, the computation of the extremal mapping 
between connected surfaces embedded in E 3 can be reduced to the computation of 
the extremal mapping between the conformal domains in C. 

In the subsequent section, we propose an algorithm, called the quasi- conformal 



(QC) iteration to solve the above minimization problems (4.5) 



5. Main Algorithm. In this section, we describe an iterative scheme, called 



the quasi- conformal (QC) iteration, for solving the variational problem (4.5). The 
QC iteration is based on the Linear Beltrami Solver(LBS). The LBS will firstly be 
explained in detail. QC iteration will then be described. 



10 



Lok Ming Lui et al. 



Practically speaking, 2D domains or surfaces in R 3 are usually represented dis- 
cretely by triangular meshes. Suppose K\ and K 2 are two surface meshes with the 
same topology representing Si and 52- We define the set of vertices on K\ and K 2 
by V 1 = {vj}f =1 and V 2 = {i> 2 }™ =1 respectively. Similarly, we define the set of tri- 
angular faces on K x and K 2 by F 1 = {Tj}JL 1 and F 2 = {T 2 }f =1 . Our goal is to 
look for a piecewise linear homeomorphism between K± and K 2 that approximates 
the Teichmuller extremal mapping between Si and S 2 . 

5.1. Linear Beltrami Solver. Our goal is to look for an optimal Beltrami co- 
efficient (BC) associated to the desired Teichmuller mapping. Every quasi-conformal 
mapping is associated to a unique BC. Given a BC, it is important to have an algo- 
rithm to reconstruct the associated quasi-conformal homeomorphism. 

Suppose / : K\ — > K 2 is an orientation preserving piecewise linear homeomor- 
phism between K\ and K 2 . We can assume K\ and K 2 are both embedded in R 2 . 
In case K\ and K 2 are surface meshes in M 3 , we first parameterize them conformally 
by 4>i : Ki -> D t C R 2 and 4> 2 : K 2 -> D 2 C R 2 . The composition of / with the 
conformal parameterizations, / := (f> 2 o f o (p^ 1 , is then an orientation preserving 
piecewise linear homeomorphism between D\ and D 2 embedded in R 2 . In this paper, 
we assume the topology of the surface mesh is either a connected open surface or a 
genus-0 closed surface. In other words, the conformal domain Dj (i = 1, 2) can either 
be a 2D rectangle, unit disk, punctual disk or unit sphere. 

To compute the quasi-conformal mapping, the key idea is to discretize Equation 
|3.7| with two linear systems. 

Given a map / = (u + \f— lu) : K\ — > K 2 , we can easily compute its associated 
Beltrami coefficient ///, which is a complex- valued function defined on each triangular 
faces of K\. To compute /ij, we simply need to approximate the partial derivatives 
at every face T. We denote them by D x f(T) = D x u + \f^\D x v and D y f{T) = 
D y u + y/—lD y v respectively. Note that / is piecewise linear. The restriction of / on 
each triangular face T can be written as: 

ti ( \ ( aTX + b T y + r T \ , , 

fMx > y)= {c TX + d T y + ST ) ^ 

Hence, D x u(T) = a T , D v u(T) = b T , D x v(T) = c T and D y v{T) = d T . Now, the 
gradient Vy/ := (D x f(T), D y f(T)Y on each face T can be computed by solving the 
linear system: 

l 1 '^ )V T fi=l f A~f,~, . (5-2) 
v 2 - v J J \ A(«HM 



where [t*o,wl] and [wo,^] are two edges on T. By solving equation 5.2 ax, &t, ct and 
dx can be obtained. The Beltrami coefficient fJ,/(T) of the triangular face T can then 
be computed from the Beltrami equation |3.1| by: 



[ax + d T ) + v— 1(ct — bx) 



Equation 3.5 and 3.6 are both satisfied on e very tri angu lar faces. Let fif(T) 

can be obtained. 



px + v — 1 t x ■ The discrete versions of Equation 3.5 and 3.6 



-d T = oti{T)ax + a 2 (T)b T ^ 
c T = a 2 {T)a T + a 3 (T)b T 
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and 



-b T = ai(T)c T + a 2 (T)d T 
a T — a 2 {T)cT + a 3 (T)d T 



where: ai (T) = a 2 (T) = -j^; a 3 (T) = ^ffi* . 

In order to discretize Equation |3.7[ we need to introduce the discrete divergence. 
The discrete divergence can be defined as follows. Let T = [v~i, Vj, vu] and wi = f(vi) 
where I — i,j or k. Suppose vj = gi + hi and wi = sj + ti (I = i,j, k). 
Using equation 5.2 et^, &t, ct and dx can be written as follows: 

a T = Ajsi + Ajs 3 + Als k - b T = Bj Sl + BJs^ + Bls k ; ^ 
c T = AfU + Aft 3 + A%t k ; d T = B?U + Bjtj + B^t k ; 

where: 

A? = (hj - h k )/Area(T), Aj = (h k - hi)/Area{T), A T k = (hi - hj)/Area{T); 

Bj = (g k - gj )/Area(T), Bj = ( 9i - g k )/Area{T), B k = ( 9j - g t )/Area(T); 

(5.7) 

Suppose V = (Vi, V 2 ) is a discrete vector field defined on every triangular faces. 
For each vertex Uj, let Ni be the collection of neighborhood faces attached to Uj. We 
define the discrete divergence Div of V as follows: 

Div(V)( Vi ) = £ AjVi{T) + BjV 2 {T) (5.8) 
TeNi 

By careful checking, one can prove that 

]T A? b T = Bfa T ; ]T Afd T = ]T Bjc T . (5.9) 

TeNi TeNi TeNi TeNi 



This gives 

Dili I w , ( , j, , t ; , ■ 

D x u I \ D x v 



Div [ ® yU \=Q and Div ( ® yV ) = (5.10) 



As a result, Equation (3.7) can be discretized: 

Div (a ( ^ xU Y) = and Div (a( „* V ) ) = U (5.11) 



D y U I I V V 



where A = [ ai a2 I . This is equivalent to: 
1 a 2 a 3 ' 



J2 Aj[ai{T)a T + a 2 (T)b T ] + Bj[a 2 {T)a T + a 3 (T)b T ] = (5.12) 



TeNi 



5^ Af[ai(T)c T + a 2 (T)d T ] + Bj[a 2 {T)c T + a 3 (T)d T ] = (5.13) 
TeNi 
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(H0| (1,0) (0,0) (1.0) 

(A) Simply-connected open surfaces 





(B) Geiius-0 closed surfaces 



Fig. 5.1. Illustration of how Beltrami representations for homeomorphisms between meshes can 
be computed. (A) shows the case of a homeomorphism between simply-connected open meshes. The 
two meshes are mapped to a unit square by harmonic parameterizations. (B) shows the case of a 
homeomorphism between genus-0 closed surface meshes. The two meshes are parameterized onto a 
triangle in R2, after cutting away a triangular face on each mesh. 



for all vertices Vi G D. Note that ax and bx can be written as a linear combination of 



the x-coordinates of the desired quasi-conformal map /. Hence, equation 5.12 gives 
us the linear systems to solve for the x-coordinate function of /. Similarly, ct and 

can also be written as a linear combination of the y-coordinates of the desired 
quasi-conformal map /. Therefore, equation |5. 13| gives us the linear systems to solve 
for the y-coordinate function of /. 

Besides, / has to satisfy certain constraints on the boundary. One common 
situation is to give the Dirichlet condition on the whole boundary. That is, for any 
v b G BK X 



f(v b ) =w b e dK 2 



(5.14) 



Note that the Dirichlet condition is not required to be enforced on the whole 
boundary. The proposed algorithm also allows free boundary condition. In the case 
that Ki and K 2 are rectangles, the desired quasi-conformal map should satisfy 



/(0) = 0;/(l) = l/(i) = i/(l + i) = l + »; 
Re(/) = on arc [0, i]; Re(/) = 1 on arc [1, 1 + 
Imag(/) = on arc [0, 1]; Imag(/) — 1 on arc [i, 1 + i] 



(5.15) 



When Ki (i — 1,2) is an unit disk, we can parameterize it onto a domain _Dj, 
which is a triangle with boundary vertices p l , p\ and p\. p\ is on the y-axis whereas p\ 
and p l 2 are on the x-axis. This can be done by removing a triangular face at the point 1 
and map Ki to the upper half plane using a Mobiiis transformation: iji(z) = \/— l-^rf • 
In this case, the desired quasi-conformal map / should satisfy 

/(pj) = pI; /(pi) = p\ and Imag(/) = on arc \pl,p{}; (5.16) 

When Ki (i = 1, 2) is a genus-0 closed surface mesh, we can again parameterize 
it onto a domain Dj, which is a triangle with boundary vertices Pq, p\ and p\. This 
can be done by removing a triangular face at the north pole and map Ki to the 2D 
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plane using stereographic projection. In this case, the desired quasi-conformal map / 
should satisfy 

f(pl) = Po; M) = Pi and f{p\) = p\ (5.17) 

Suppose landmark correspondences {pi}" =1 ■<-> {<Zi}" =1 are enforced, one should 
add this constraint to the linear system. Mathematically, it is described as f(pi) = Qi 
(z = l,2,...,n). 

Equations |5 . 1 2| and |5 . 13| together with the above boundary conditions give a non- 
singular linear system to solve for /. The linear system is symmetric positive definite. 
Hence, it can be solved effectively by the conjugate gradient method. We call this 
algorithm the Linear Beltrami 5o/uer(LBS). Given a Beltrami coefficient v, we denote 
the obtained quasi-conformal map from LBS by LBS(f). If landmark constraints are 
enforced, we denote it by LBS/,m(z/). 

We note that given an arbitrary Beltrami coefficient v and arbitrary landmark 
correspondences, a quasi-conformal mapping associated to v might not exist. How- 
ever, the Linear Beltrami Solver looks for the best quasi-conformal mapping whose 
Beltrami coefficient closely resemble to v. 

5.2. Quasi-conformal(QC) iterations. With the Linear Beltrami Solver, one 
can easily obtain the best quasi-conformal mapping associated with a given BC. In 
order to obtain the extremal mapping /, our goal is to iteratively search for the 
optimal BC associated to /. With the optimal BC, the desired extremal mapping / 
can be easily reconstructed using the Linear Beltrami Solver. 

Recall that our problem of computing the extremal mapping can be converted 
into an optimization problem: 

(u, f) = argmin^^JI \v\ U + 1 1 W| | 2 } (5.18) 

subject to: (1) v — fi(f) and IMloo < 1; (2) v = k — for some constant k and 
holomorphic function ip : D\ — > C; and (3) / satisfies certain boundary condition 
and/or landmark constraints. Note that the boundary condition in (3) can either be a 
Dirichlct condition defined on the whole boundary or free boundary condition. In this 
subsection, we introduce the Quasi- conformalf QC) iteration to solve the optimization 
problem. 

The QC iteration starts with an initial map f : D\ — > D 2 satisfying the given 
boundary condition and landmark constraints. The initial map is chosen to be the 
quasi-conformal mapping obtained from LBS associated to the BC fio — 0. In other 
words, 

fo = LBS iM (Mo := 0) (5.19) 

Note that with the enforced landmark constraints, the Beltrami coefficient asso- 
ciated to fo might not be equal to i>q. The Linear Beltrami Solver simply look for 
the best quasi-conformal mapping whose Beltrami coefficient resemble to v as much 
as possible. Let Uq be the Beltrami coefficient associated to fo. This gives us a pair 
(fo, vo), for which v = fi(fo)- 

Now, in order to minimize the energy function E 2 satisfying condition (2), we 
propose to perform a Laplace smooth £ and averaging 21 on vq- The Laplace smooth 
£, which aims to minimize E 2 , is given by the following: 



Z(v Q )(T) := ^ /|Nbhd(T)| 

T ; eNbhd(T) 



(5.20) 
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where T is a triangular face of K\, Nbhd(T) is the set of neighborhood faces of T 
and |Nbhd(T)| is the number of neighborhood faces in the set Nbhd(T). Set li\{T) = 



A aims to obtain an optimal v satisfying the condition (3) in the optimization problem. 
An updated quasi-conformal, f\ can then be obtained by LBS: f\ = LBSlm(ah)- And 
an updated Beltrami coefficient, v\ := can be computed. Thus, we get a new 

pair (fi,vi). 

The procedure continues until the iteration converges. More specifically, given 
the pair (f n ,u n ) obtained at the n iteration, we can obtain a new pair (f n +i, v n +i) 
as follows: 



Consequently we get a sequence of pair (/„, u n ), which converges to the optimal 
Beltrami coefficient associated to the extremal mapping or optimized Teichmuller 
mapping. In practice, we stop the iteration when \\v n +i ~ u n\\ < e - 

We summarize the QC iteration as follows. 

Algorithm 5.1 : (QC iteration for open surfaces) 

Input : Triangular meshes: K\ and K 2 and the desired boundary condition 
Output : Optimal Beltrami coefficient v and the extremal mapping f 

1. Obtain the initial mapping f = LBSlm(mo : = 0). Set v a = /li(/o); 

2. Given v n , compute /U n +i := A(£(v n )); Compute := LBS^mIMk+i) an d 
set v n+1 := (j,(f n+ i); 

3. 7/||f„ + i — v n \\ > e, continue. Otherwise, stop the iteration. 

The QC iteration can also be applied to the case when Z)j (i = 1, 2) is a unit 
sphere. In other words, given a set of landmark constraints between the unit sphere, 
our goal is to look for the Teichmuller extremal mapping / : D\ —¥ D%. However, 
special attention has to be paid in this case. 

Denote the landmark correspondence by {Pi}™ =1 -H- {<7i}" =1 . We can assume that 
the north pole is fixed. If not, it can also be achieved by a Mobius transformation. 
The LBS can be applied to unit spheres, by stereographically projecting Di onto a 
big triangles in M 2 . However, numerical error near the north pole is inevitable. We 
therefore propose an alternating scheme to fix this problem. 

For the initial map, we add the vertices near the north pole {rij}J l =1 (z > 0.99) as 
landmarks and fix all {nj}JL 1 . We then compute the Teichmuller mapping f using 
Algorithm 5.1. Numerical error will be introduced near the north pole. To fix it, in 
our next step, we consider the vertices {sj}JL 1 near the south pole (z < —0.99) as 
landmarks. The correspondence is given by: Sj o fi(sj). Rotate the south pole of D, 
to the north pole by a Mobius transformation. We can again compute the Teichmuller 
mapping f\ using Algorithm 5.1. 

We continue this process until the iteration converges. More specifically, at the 
n iteration where n is an even integer, we add vertices {sj}"L 1 around south pole 




(HiT) = A(jl!){T) :=( 



all faces of K\ 

No. of faces of K x V~i( T )l 



(5.21) 



H n+1 := A{£{v n )); 
fn+i '■= LBS LM (/i„ +1 ); 

V n +1 ■= Kfn+l)- 



(5.22) 
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Fig. 5.2. Example of the Teichmiiller extremal mapping of the disk with fixed Dirichlet boundary 
condition. (A) shows the boundary condition. (B) shows the energy and the sup-norm of the BC in 
the QC iteration. (C) shows the obtained extremal mapping, visualized using the texture mapping. 
(D) shows the histogram of the BC norm. 



as landmarks. Set correspondence as: Sj ■<-» f n (sj). Rotate the south pole of Di to 
the north pole by a Mobiiis transformation, and obtain the Teichmiiller mapping / n +i 
using Algorithm 5.1. When n is an odd integer, we add vertices {n-j}™ =1 around north 
pole as landmarks. Set correspondence as: rij <-$■ f n {nj) and obtain the Teichmiiller 
mapping f n+1 using Algorithm 5.1. Set v n+ x = fi(f n )- 

This alternating process between the north pole and the south pole continues 
until H^n+i -Vn\\ < e< 

When Di (i = 1,2) is a unit disk, the LBS would also introduce numerical error 
near 1. To fix it, the same alternating algorithm between 1 and -1 can be applied. 

The detailed algorithm can be summarized as follows: 



Algorithm 5.2 : (QC iteration for genus-0 closed surfaces) 

Input : Triangular meshes: K\ and K2 and the desired boundary condition 

Output : Optimal Beltrami coefficient v and the extremal mapping f 

1. Add vertices around north pole as landmarks and fix their positions. Obtain 
the initial Teichmiiller mapping fo using Algorithm 5.1. Set v$ = /x(/o); 

2. Given f n and v n . When n is even, add vertices {sj around south pole 
as landmarks. Set correspondence as: Sj f-> f n (sj)- Rotate the south pole of 
Di to the north pole. When n is odd, add vertices {sj}™ =1 around south pole 
as landmarks. Set correspondence as: Sj O f n (sj). Obtain the Teichmiiller 
mapping f n +i using Algorithm 5.1. Set v n+ \ = /i(/ n ) ; '; 

3- If \\v n +\ — v n \ \ > e, continue. Otherwise, stop the iteration. 



6. Numerical experiments. In this section, we evaluate our proposed algo- 
rithm numerically by synthetic examples. 
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Fig. 5.3. Another example of the Teichmuller extremal mapping of the disk with fixed Dirichlet 
boundary condition. (A) shows the obtained extremal mapping, visualized using the texture mapping. 
(B) shows the histogram of the 
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Fig. 5.4. Teichmuller extremal mapping of the disk with only 8 landmark points constraints 
on the boundary. (A) shows the Teichmuller extremal mapping. (B) shows the histogram of the 
BC norm. (C) shows the energy and the sup-norm of the BC in the QC iteration. (D) shows the 
automatically obtained optimal boundary correspondence. (E) shows the histogram of the BC norm 
under harmonic map with arc-length parameterized boundary condition. (F) shows the histogram of 
the BC norm under the Teichmuller mapping with arc-length parameterized boundary condition. 



6.1. Extremal mapping of simply-connected domains. In our first numer- 
ical experiment, we test our method to compute the extremal mapping of the disk 
with a given Dirichlet boundary condition. A Dirchlet condition on the whole bound- 
ary is given as shown in Figure 5.2 'A). Using the QC iteration, we iteratively obtain 



the pair (f n ,y„ = As shown in Figure |5^2[B), E 2 {f n ,v n ), \Wn\\oo and ||Vi/|| 2 

all decrease as iteration increases. The resulting Teichmuller extremal mapping is as 
shown in Figure [5~2^ C), which is visualized using the texture mapping. Note that the 
original texture is deformed under the extremal mapping. However, the dilations of 
the ellipses deformed from the small circles are the same. It means the norm of the 
BC is constant everywhere. The histogram of the norm of the BC is also shown in 
Figure 5.2 D), which again demonstrates the norm of the BC is equal to a constant 
k = 0.15016. The standard deviation of the BC norm is 0.0034373. Note that the 
Dirichlet boundary condition can be of arbitrary shapes. Figure [573] shows the result 
of finding the Teichmuller extremal mapping between the disk and the ameba with 



Teichmiiller mapping and applications 



17 




(A) 



CB) 



4?o*o°o*c£c 
b # o°o # o°o*d 



(C) 



it- 



(D) 



Fig. 6.1. Teichmiiller extremal mapping of the annulus with fixed Dirichlet boundary condition. 
(A) shows the boundary condition. (B) shows the energy and the sup-norm of the BC in the QC 
iteration. (C) shows the obtained extremal mapping, visualized using the texture mapping. (D) 
shows the histogram of the BC norm. 




Fig. 6.2. Teichmiiller extremal mapping of the multiply- connected domain containing three holes 
with fixed Dirichlet boundary condition. (A) shows the Teichmiiller mapping Dirichlet boundary 
condition. (B) shows the histogram of the norm of the BC. (C) shows the energy, the sup-norm of 
the BC and the L2 norm of the gradient of the nom of BC in the QC iteration. 



given boundary conditions. The resultant norm of BC is also constant everywhere, as 
demonstrated in Figure [53] (B). 

Our algorithm also applies to situation when only a few landmark constraints are 
enforced on the boundary (instead of the Dirichlet condition defined on the whole 
boundary) . In Figure |5.4| we test our algorithm to compute Teichmiiller extremal 
mapping of the disk with only 8 landmark points constraints on the boundary. (A) 
shows the the Teichmiiller extremal mapping. Again, the dilations of the ellipses 
deformed from the small circles are the same, meaning that the norm of the BC 
is constant everywhere. (B) shows the histogram of the BC norm. The norm k 
of the BC is equal to 0.201. (C) shows the energy and the sup-norm of the BC 
in the QC iteration. Our algorithm also automatically detect the optimal boundary 
correspondence. (D) shows the obtained optimal boundary correspondence. (E) shows 
the histogram of the BC norm under harmonic map with arc-length correspondence 
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Fig. 6.3. Teichmiiller extremal mapping of the multiply- connected domain containing six holes 
with fixed Dirichlet boundary condition. (A) shows the Teichmiiller mapping Dirichlet boundary 
condition. (B) shows the histogram of the norm of the BC. 



on the boundary. Note that the distribution of the conformality distortion is highly 
non-uniform. (F) shows the histogram of the BC norm under the Teichmiiller mapping 
with arc- length parameterized boundary condition (of which the landmark constraints 
are satisfied). Although a Teichmiiller mapping can still be obtained, the norm of the 
BC is equal to 0.23 which is higher than the case when only 8 points landmark 
constraints are enforced. Hence, the obtained Teichmiiller mapping is not extremal. 

6.2. Extremal mapping of multiply-connected domains. Our method can 
also be applied to multiply-connected domain. In Figure |6.1[ we test our method 
to compute the extremal mapping of an annulus with Dirichlet boundary condition. 
The boundary condition is given in (B). Again, the energy functional Ei is decreasing 
under the QC iteration. The extremal mapping is as shown in (C). The mapping is a 
Teichmiiller map, since the BC norm is constant everywhere as shown in (D). 



Figure 6.2 and 6.3 show the results of finding the Teichmiiller extremal map- 
ping with 3 and 6 holes respectively. Figure |6.2| shows the Teichmiiller map. The 
corresponding BC norm, sup-norm of BC and the L2 norm of the gradient of the 



norm of BC are shown in (B) and (C) respectively. Figure 6.3 demonstrates the re- 
sult of a more complicated multiply-connected domain having 6 holes. These results 
show that our algorithm can again compute the Teichmiiller Extremal mapping of 
multiply-connected domains efficiently and accurately. 

6.3. Optimized Teichmiiller mapping with interior landmark constraints. 

Our algorithm can compute an optimized Teichmiiller mapping with interior landmark 
constraints enforced. Figure |6.4| shows the optimized Teichmiiller mapping between 
the disk with 25 interior landmark constraints enforced. (A) shows the 24 landmark 
constraints. (B) shows the obtained Teichmiiller mapping, visualized using the tex- 
ture mapping. (C) shows the energy, the sup-norm of the BC and the L-2 norm of the 
gradient of the BC in the QC iteration. Note that the sup-norm of the BC decreases 
during the iteration, which illustrates the mapping converges to an optimal mapping 
minimizing the conformality distortion. (D) shows the histogram of the norm of the 
BC. The norm k of the BC is uniformly equal to 0.2. (E) shows the histogram of the 
norm of BC with arc-length boundary correspondence enforced. The norm k of the 
BC is equal to 0.28, which means the Teichmiiller mapping is not an extremal one. 

We have test our algorithm to compute the optimized Teichmiiller mapping of the 
disk with 3 interior landmark curves constraints enforced in Figure |6.5| (A) shows 
the Teichmiiller mapping with 3 landmark curves constraints enforced. (B) shows the 
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Fig. 6.4. Teichmiiller mapping between the disks with 24 interior landmark constraints enforced. 
(A) shows the 24 landmark constraints. (B) shows the Teichmiiller mapping. (C) shows the energy, 
the sup-norm of the BC and the L-2 norm of the gradient of the BC in the QC iteration. (D) 
shows the histogram of the norm of the BC. (E) shows the histogram of the BC norm under the 
Teichmiiller mapping with arc-length parameterized boundary condition. 



energy, the sup-norm of the BC and the L-2 norm of the gradient of the BC in the 
QC iteration. The energy is decreasing, indicating that the algorithm converges to 
an optimized Teichmiiller mapping. (C) shows the histogram of the norm of the BC. 
The norm of the BC is accumulated at 0.53. 

We also test our algorithm to the case when we only have the correspondence of 
the interior landmarks. Figure [6~6| shows an example of finding Teichmiiller extremal 
mapping with 20 interior landmark points only. (A) shows the constraint of the 
feature points. (B) shows the Teichmiiller mapping. (C) shows the histogram of the 
norm of the BC. It shows that even boundary condition is not provided, our algorithm 
can still obtain both the Teichmiiller extremal mapping which satisfies the landmark 
constraints and the corresponding optimal boundary condition. 

Our algorithm also applies to computing the optimized Teichmiiller extremal map- 



ping between the unit sphere with interior landmark constraints enforced. Figure 6.7 
shows an example of Teichmiiller extremal mapping with 10 interior landmark con- 
straints enforced. (A) shows landmark constraints on the sphere. (B) shows the 
Teichmiiller mapping. (C) shows the energy, the sup-norm of the BC and the L-2 
norm of the gradient of the BC in the QC iteration. Again, the energy is decreasing, 
indicating that the algorithm converges to an optimized Teichmiiller mapping. (D) 
shows the histogram of the norm of the BC. The norm of the BC is concentrated near 
0.21, meaning that the mapping is a Teichmiiller mapping. 



7. Applications. In this section, we apply our proposed algorithms for com- 
puting landmark matching Teichmiiller mappings to practical problems. More specif- 
ically, we will consider the problems of computing brain landmark matching registra- 
tions, constrained texture mappings and human face registrations. 
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Fig. 6.5. Teichmiiller mapping between the disks with 3 interior landmark curves constraints 
enforced. (A) shows the Teichmiiller mapping with 3 landmark curves constraints enforced. (B) 
shows the energy, the sup-norm of the BC and the L-2 norm of the gradient of the BC in the QC 
iteration. ( C) shows the histogram of the norm of the BC. 




Fig. 6.6. Teichmiiller mapping between the disks with 20 interior landmark constraints en- 
forced. (A) shows the 20 landmark constraints. (B) shows the Teichmiiller mapping. (C) shows the 
histogram of the norm of the BC. 



7.1. Brain landmark matching registration. Landmark-based surface reg- 
istrations are commonly applied for finding meaningful 1-1 correspondences between 
human brain cortical surfaces. On cortical surfaces, sulcal landmarks can be labeled 
either manually by neuroscientists or automatically based on various geometric quanti- 
ties. The sulcal landmarks are important anatomical features. It is therefore desirable 
to obtain a registration between the cortical surfaces with least geometric distortion, 
which matches the sulcal landmarks as much as possible. Our algorithms for com- 
puting landmark matching Teichmiiller mappings can be applied. In Figure |7.1[ we 
apply our algorithm to compute the Teichmiiller mapping between 2 different brain 
surfaces with 3 corresponding landmarks labeled. (A) shows the corresponding sulcal 
landmarks, indicated by different colors. (B) shows the obtained Teichmiiller extremal 
mapping with 3 landmark constraints enforced, visualized by the circle packing tex- 
tures. The sulcal landmarks are exactly matched under the mapping. (C) shows the 
histogram of the norm of the associated BC. The norm is a constant showing that the 
obtained registration is indeed a Teichmiiller mapping. We also test the method to 
register cortical surfaces with more sulcal landmarks. In Figure |7.2[ we compute the 
Teichmiiller extremal mapping between 2 brain surfaces with 6 corresponding sulcal 
landmarks labeled. The obtained registration and the norm of its associated BC is 
shown in (B) and (C) respectively. The landmarks are exactly matched. Again, the 
norm of the BC is a constant, showing that the obtained registration is a Teichmiiller 
mapping which minimizes the conformality distortion. 

7.2. Constrained texture mapping. Texture mapping is one of the major 
photorealistic techniques in computer graphics to generate realistic and visually rich 
3D surfaces 23 , 24 . It is usually done by putting each surface mesh in correspondence 



Teichmiiller mapping and applications 



21 



(A) (B) 



"a j\ 



(Q (D) 

Fig. 6.7. Teichmiiller mapping between the spheres. (A) shows landmark constraints on the 
sphere. (B) shows the Teichmiiller mapping. (C) shows the energy, the sup-norm of the BC and 
the L-2 norm of the gradient of the BC in the QC iteration. (D) shows the histogram of the norm 
of the BC. 

with a 2D image [2j [30l EH ESI EB ES [25] . Such a correspondence between the surface 
mesh and the image is called the texture mapping. Constrained texture mappings 
are popularly used, in which the texture mappings are guided by landmark features 
labeled interactively by users. Ideally, the texture mapping should have minimum 
distortion, while matching the landmark points exactly. We apply our algorithms to 
compute the landmark matching Teichmiiller extremal mapping between the surface 
mesh and the image, and use it as the texture mapping. In Figure |7.3[ we map a cat 
image onto the human face surface. (A) shows the corresponding landmark points 
labeled manually on the human face and the texture image. The landmark matching 
Teichmiiller extremal mapping is computed, and is used as texture mapping to project 
the cat image onto the human face surface. The textured surface is as shown in 
(B). (C) shows the norm of the associated BC of the texture mapping. The norm 
is approximately a constant. It means the texture mapping computed is indeed a 
Teichmiiller extremal mapping, which minimizes the conformality distortion. 

In Figure [7^4} we further test our algorithm on a multiply-connected human face. 
(A) and (B) shows the texture(tiger) image and a multiply-connected human face. 
Corresponding landmark points are labeled manually on the texture image and the 
surface mesh. In (C), the surface mesh is mapped to a multiply-connected domain in 
2D by a Teichmiiller extremal mapping matching the landmark points exactly. (D) 
shows the textured surface. (E) shows the norm of the associated BC of the texture 
mapping. The norm is approximately a constant, which means the texture mapping 
computed is indeed a Teichmiiller extremal mapping. 

7.3. Human face registration. In face recognition, finding accurate spatial 
correspondences between human faces is an a crucial process to compare and recog- 
nize faces effectively. Corresponding features can be extracted on human face based on 



22 



Lok Ming Lui et al. 




Fig. 7.1. (A) shows 2 brain surfaces with 3 corresponding landmarks. (B) shows the Te- 
ichmiiller extremal mapping with 3 landmark constraints enforced. (C) shows the histogram of the 
norm of BC. 

curvatures, such as high curvature points near nose tips and lips. Accurate face regis- 
tration can then be obtained by computing a mapping that matches the corresponding 
features. Landmark matching Teichmiiller extremal mapping, which minimizes the 
geometric distortion, can then be used. In Figure |7.5[ we apply our algorithm to 
compute the registration between a male and female human faces. The human faces 
are both simply-connected open surfaces. Corresponding feature points are labeled 
on both faces. The obtained Teichmiiller mapping is obtained, which is visualized by 
texture mapping. The corresponding features are exactly matched. (C) shows the 
histogram of the norm of the BC, which is almost a constant. This demonstrates the 
obtained registration is a Teichmiiller mapping. 

Our algorithm can also be applied to obtain registration between multiply-connected 
human faces. Figure |7l)| shows two multiply-connected human faces. Corresponding 
feature landmarks are labeled. Teichmiiller extremal mapping matching the features 
exactly is computed, as shown in (B). It is again visualized by texture mapping. (C) 
shows the histogram of the norm of the BC. Again, it is almost a constant, which 
demonstrates that obtained registration is a Teichmullcr mapping minimizing the 
conformality distortion. 

8. Conclusion. We address the problem of computing Teichmiiller extremal 
mapping between surfaces, which minimizes the maximal conformality distortion. 
The proposed algorithm can be applied to obtain a landmark matching registration 
between surface meshes. Given a set of corresponding landmark points or curves 
defined on both surfaces, a unique landmark matching quasi-conformal registration 
can be obtained, which minimizes the conformality distortion. In this paper, we 
propose an efficient iterative algorithm, called the Quasi-conformal (QC) iterations, 
to compute the Teichmiiller extremal mapping. The key idea is to represent the set 
of diffeomorphisms by Beltrami coefficients (BCs). We then look for an optimal BC 
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Fig. 7.3. (A) shows a human face and a texture image of a cat. Corresponding landmark points 
are labeled on the surface and the texture image. We compute the Teichmiiller extremal mapping 
that matches the landmark points. The Teichmiiller mapping is used as constrained texture mapping 
to project the texture image onto the surface, as shown in (B). (C) shows the histogram of the norm 
ofBC. 



associated to the desired Teichmiiller mapping. The associated diffeomorphism can be 
efficiently reconstructed from the optimal BC using the Linear Beltrami Solver(LBS). 
Using our proposed method, the Teichmiiller mapping can be accurately and efficiently 
computed within 10 seconds. The obtained registration is guaranteed to be bijective. 
Besides, Teichmiiller mapping with soft landmark constraints can also be computed 
using our proposed algorithm. It becomes useful when landmark features cannot be 
accurately located, and hence it is better to compute registration with landmarks 
approximately (but not exactly) matched. We applied the proposed algorithm to 
real applications, such as brain landmark matching registration, constrained texture 
mapping and human face registration. Experimental results shows that our method 
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Fig. 7.4. (A) a texture image of a tiger. (B) shows a multiply-connected human face. Cor- 
responding landmark points are labeled on the surface and the texture image. In (C), the surface 
mesh is mapped to a multiply- connected domain in 2D by a Teichmiiller extremal mapping matching 
the landmark points exactly. The Teichmiiller mapping is used as constrained texture mapping to 
project the texture image onto the surface, as shown in (D). (E) shows the histogram of the norm 
ofBC. 






Fig. 7.5. Teichmiiller extremal mapping of the simply- connected domain with landmark point 
constraints (A) shows the two faces with landmark point constraints. (B) shows the Teichmiiller 
extremal mapping of the two faces. The resultant mapping is illustrated by texture mapping. (C) 
shows the histogram of the norm of BC. 



is effective in computing a non-overlap landmark matching registration with least 
amount of conformality distortion. 
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